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Abstract 

o : 

In this paper we study the structure of the manifold of solitary waves in some deformations of 50(2) 
symmetric two-component scalar field theoretical models in two-dimensional Minkowski space. The 
deformation is chosen in order to make the analogous mechanical system Hamilton- Jacobi separable 
in polar coordinates and displays a singularity at the origin of the internal plane. The existence of the 
singularity confers interesting and intriguing properties to the solitary waves or kink solutions. 

in 

C ' 1 Introduction 



> 



Solitary waves are at the heart of astonishing phenomena in diverse physical systems that can be modeled 
by non-linear equations. For instance, they describe the behavior of interfaces in magnetic materials pQ 
and in ferroelectric crystals [2] in Condensed Matter. Solitary waves have biotechnical and biomedical 



applications [3] and also seed the formation of structures in Cosmology [4j. In several disguises, kinks, 
topological defects, domain walls, membranes, solitary waves arise in many branches of Theoretical 
Physics [21 E]- For this reason, the search for solitary waves in some types of PDE is an active topic 
in non-linear science. Among these equations we find the non-linear Klein-Gordon equation, extensively 
cited in the physical literature 0[S]. This equation, generalized to iV-real scalar fields, reads as follows: 

d 2 <j> a d 2 <P a ,dU 

where U(4>i, . . . , 0tv) is a potential function of the scalar fields <j> a , whereas Jjjp are non-linear functions 
^ ■ of the fields. PDE (JT]) can be understood as the Euler-Lagrange equations associated with the functional 



action 



S = I d 2 x 



1 N 

-Y^d^ a d»4> a -U{^...,cp N ) 

0=1 



governing the dynamics of a (1+1) dimensional scalar field theory. In this framework, a solitary wave 
is a localized non-singular solution of the non-linear field equation HP whose energy density, as well as 
being localized, has space-time dependence of the form: e(t,x) = e(x — vt), where v is some velocity vector 
according to Rajaraman [8]. Use of Lorentz invariance allows us to investigate the existence of some 
kinds of solitary waves by reducing PDE ([T]) to the following ODE: 
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Therefore, the search for solitary waves or kinks, finite energy solutions to the static field equations ([2]), is 
tantamount to solving a analogous mechanical problem for a unit-mass point particle moving in a plane 
with coordinates <j) = {fa, fa) under the influence of a potential —U(4>), if the variable x plays the role 
of "time" 0. 

The sine-Gordon and 4 models, profusely dealt with in the literature, are the basic examples of (TjQ), 
respectively governed by the PDE equations: 

d 2 4> d 2 d) d 2 4> d 2 4> , 2 , 



dt 2 dx 2 dt 2 dx 2 

The corresponding potential terms are U(4>) = 1 — cos0 and U{cf)) = \(<j) 2 — l) 2 and for both systems 
solitary or traveling wave solutions exist: the well known sine-Gordon soliton and A(c/>)| kink. The 
peculiar non-dispersive character of these non-linear waves is related to the structure of the set Ai of 
zeroes of the potential U(fa. Ai is a discrete set with more than one element. The field profile of these 
solitary waves connects two elements of Ai asymptotically. For instance, the kink 4>k( x ) = ^tanhx in the 
cj) model interpolates between the two zeroes cj>± = ±1 of U(fa: 4>k(x = —00) = ±1, 4>k{ x = 00) = =Fl. 
Obviously, the greater the number of elements in Ai, the richer the solitary wave variety. Thus, other 
models have been considered in the literature in the search for these kinds of non-linear waves, such as 
the (p 6 model with potential U(fa = \4> 2 {4> 2 - l) 2 , the 8 model with U(<f>) = \{4> 2 - l) 2 {4> 2 - a 2 ) 2 in 
the one-component scalar field theory framework [9] and, more recently, polynomial interactions of any 
order built from Chebyshev polynomials |lUj . 

Another way of obtaining a richer structure in the set Ai, and hence in the solitary wave manifold, 
is to increase the dimension of the internal space, i.e, by considering theories with more scalar field 
components. This is an important qualitative step, as noted by Rajaraman [8]: This already brings us to 
the stage where no general methods are available for obtaining all localized static solutions, given the field 
equations. However, some solutions, but by no means all, can be obtained for a class of such Lagrangians 
using a little trial and error. 

Straightforward generalization of the c/> 4 and </> 8 models to two-component scalar field theory leads to 
the potential energy densities: 

v^(fa,fa) = (rt + rt-l) 2 (3) 
v^{fa,fa) = (0 2 + 2 -l) 2 (0 2 + ^-a 2 ) 2 . (4) 

The action functional is invariant under SO{2) rotations in the <p\ — 4>2 internal plane. The zeroes of 
v^ A ' and y( B \ however, are not invariant under the SO(2) action and the orbits, Ai, are continuous 
manifolds in these cases: 

M (A) = S R=i^ M (B) = s r=i u S R=a^ respective i y . Goldstone b osons arise 
in the process of quantification in this situation. Coleman proved in [Tl] that there are no Goldstone 
bosons in a sensible scalar field theory on the line. The infrared asymptotic behavior of quantum theory 
would require the modification of the potentials and 

V (B) i n 

such a way that their manifold of zeroes 
becomes a discrete set. Thus, the effect of quantum fluctuations is to add perturbations to the potentials 
© and 

U(<h,<h) = v^ B \4>iA2) + w(<t>i,<j>2) 

such that the SO (2) symmetry is explicitly broken down into discrete subgroups -acting on the new 
discrete sets of zeroes- whereas the U((f>i, fa) potential energy density remains a non-negative expression. 

There are some models in the literature that match these features. The best studied is the MSTB 
model, with potential energy density: 

£4istb (01,^2) = v^ A) (<j)i,fa) +a 2 4> 2 2 , 

where ex is a non-dimensional parameter. It was first proposed by Montonen [T2] and Sarker, Trullinger 
and Bishop [13]. Rajaraman and Weinberg pL4j discovered two kinds of kinks by the trial orbit method; 
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the TK1 (one-topological kinks) - tracing a straight line trajectory in the analogous mechanical system- 
and the TK2 (two-component topological kinks) - running through semi-elliptic orbits in the mechanical 
system. These enquiries were followed by numerical analysis |15[ [To] and with this method Subbaswamy 
and Trullinger found the existence of a whole family of two-component non-top ological kinks, which 
were named as NTK. Moreover, they discovered an unexpected fact: the NTK energy is equal to the 
addition of the TK1 and TK2 energies. This relation is known as the kink mass "sum rule" . In 1985, 
Ito explained all of these issues analytically. The crux of the matter is the separability of the Hamilton- 
Jacobi equation using elliptic coordinates in the analogous mechanical system |17[ [THj [T9"l [20] . Unlike 
models with only one scalar field, systems with two or more scalar fields have analogous mechanical 
systems that are generically non-integrable. The analogous mechanical system of the MSTB model is a 
completely integrable Type I Liouville model - separable in elliptic coordinates- and all the solitary waves 
can be found analytically . Other models exhibiting similar properties have been addressed in References 
|21 [ 1221 [2"5] . The generalization of this kind of model to three-component scalar field theory is studied in 



The goal of this work is to identify the variety of solitary wave or kink solutions in a broad family of 
models arising from perturbations of the SO (2) symmetric two-component scalar field t/> 4 and models 
([3]) and (j3|). The strategy will be to deal with analogous mechanical systems of Liouville Type II, i. e., 
with Hamilton-Jacobi equations separable in polar coordinates. These deformations necessarily involve 
a singularity at the origin of the configuration space in the mechanical system, equivalently, at the origin 
of the internal plane in the field theoretical model. Strictly speaking, the perturbation does not exist 
as a proper function when (f>± = (f>2 = because the limit of w((fii,4>2) when both <p\ and 4>2 are zero is 
either or oo, depending on the path followed to reach u>(0,0). This mathematical pathology confers 
new and intriguing properties to solitary wave solutions. For example, certain kink profiles connecting 
identical vacuum points in internal space asymptotically cannot be deformed into each other even though 
they belong to the same topological sector in the configuration space. Trajectories moving away from the 
origin involve infinite energy and some kink solutions behave as strings pinned at one point in internal 
space. Nevertheless, we shall show that in the four models that we discussed here there is a plethora 
of kink or solitary wave solutions, including topological and non-topological kinks, one-parametric kink 
families, and singular solutions. 

The organization of the paper is as follows: In Section §2 we shall describe the models to be studied 
from a generic point of view. We shall discuss the general properties of these models and state three 
interesting results. In Sections §3, §4, §5 and §6 we shall deal with the particular models that correspond 
to two perturbations of ([3]) and Applying the procedure explained in Section §2, we shall describe 
in detail the solitary wave solutions arising in these systems and unveil their hidden structure. 

2 Generalities 

In the models that we shall study, the (non-dimensional) scalar fields, 

(j}(x ,xi) = (4> 1 (x ,xi),(j)2(xo,xi)) : R 1 ' 1 -> R 2 

are maps from (l+l)-dimensional Minkowskian space-time R 1 ' to the R 2 internal space. The action 
functional 



is invariant under the Poincare transformations acting on R 1,1 , whereas the remaining symmetry trans- 
formations of our models belong to the subgroup of SO (2) rotations in the internal plane R 2 that do 
not change U(4>i,cj)2)- Our convention for the metric tensor components in Minkowski space R 1 ' 1 is 



[21 us m\. 




(5) 
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goo = —gix = 1, 512 = 521 = and only non-dimensional parameters will be considered throughout the 
paper. 

A "point" in the configuration space of the system is a configuration of the field of finite energy; i.e., 
a picture of the field at a fixed time such that the energy E, the integral over the real line of the energy 
density, 

1 fd(t)i\ 2 1 /# 2 X 2 



oo 



E W\ = / dx£[<P] , = -Jj±\ +U{<h,<h) (6) 



is finite. Thus, the configuration space is the set of continuous maps from R to R 2 of finite energy: 

C = <|V0) 6 Maps(R,R 2 )/£ < 00} 
In order to belong to C, each configuration must comply with the asymptotic conditions 

lim e M lim ^ = (7) 

x— >±oo x^±oo dx 

where Ai is the set of zeroes (minima) of the potential term t/(0). 

We shall analyze four different models, deforming the 4 and 8 potential energy densities ^ and 
dH) by two classes of perturbations: 

r 2 



where a is a non-dimensional coupling constant of the system. Thus, we consider four distinct potential 
energy densities ?7(0i,02) labelled as follows: 

C/ (/J) (01,0 2 ) =V {I) (01,02) W J) (01,0 2 ) , (10) 

where I = A,B and J = 1, 2. 

An important remark should be made here: both and u/ 2 -* are singular at the origin. The limit of 
w/ 2 )(0i, 02) when 0i — > 0, 02 — > is always infinite, independently of the path chosen in R 2 to approach 
the origin, lim^ ^ 2 )-»(o,o) w {^li 02) is, however, finite, zero, if the origin is approached through the 
abscissa axis 02 = 0, but it is equal to infinity for any other approaching path. Therefore, field profiles 
passing through the origin have infinite energy and are excluded from the configuration space in all of 
four cases, except field configurations such that all their derivatives along the 02 axis at the origin vanish 
when the perturbation chosen is . 

The invariance of the system of ODE ([2]) under spatial translations, x — > x — xq and reflections 
x — > —x makes it convenient to abbreviate the notation: x = (— l) s (x — xq), xq S R, 5 = 0, 1. In this 
manner, we shall describe a whole family of kinks and anti-kinks with their centers located at arbitrary 
points on the line, because these symmetry transformations bring solutions into solutions. 

The use of polar coordinates in the internal space R 2 is suggested by the SO (2) invariance of 
and v( B ^ as well as by the choice of perturbations and w^ 2 \ This system of coordinates is defined 
by the diffeomorphism 

(: 1+xS 1 — ► R 2 f C*(0i) = R cos if 

(Ro,<p ) — » (0?,0°) i C(<h) = Rsm<p 



< R = +v/0i + 0| < 00 > < <P = arctan^ 2 - < 2vr 
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By using polar coordinates, v^ A \ v^ B \ w^ 1 ' and become: 

f^(R) = Cv {A) = (R 2 - l) 2 ; /^(R) = Cv {B) = (R 2 - l) 2 (R 2 - a) 2 ; 

1 (1)/ \ >* (1) ° 2 ■ 2 1 (2)/ \ a* (2) ° 2 ■ 2 ¥ 

-^9 [ >W) = C w { ' = ^ sm f 5 -jp9 W) = C ™ l j = ^2 sm 2 

Thus, the first summand in (I10p is a function f^(R) that depends only on the radial variable. The 
second summand in (jlOp is the product of a function g^((p) depending only on the angular variable 
times the square of the inverse of R: 

r([7 (/j) ) = / (/) (jR) + _L 5 (j) M . 

We shall refer to this class of systems as Type II Liouville models after the type of their analogous 
integrable mechanical systems |27| . We see the reason for the inevitability of the singularity at R = 0: 
the factor in the second summand of U^ IJ \R, <p). In , however, the singularity is not seen if the 
origin is reached through the paths (p = tt and (p = 0. In the other case, w^ 2 \ only passing through the 
origin via tp = kills the singularity, but there is no continuous way out and the "particle" would become 
entrapped by the singularity. 

The general structure of the set A4 of zeroes of the potential energy density is easy to unveil: C*(U) 
vanishes if and only if the two summands on the right hand side of formula (jTTJ) are zero. Thus, 

CM = {(Ri,(fj) G R + x S 1 / f^(R%) = , g( J \tpj) = , i = 1, ...,M , j = I, N} 

and the M x N points (Ri,(fj) G Q*M lie on the knots of a lattice where the laths are the two sets of 
perpendicular straight lines: 

r R = {R = Ri,i = l,...,M} and r v = {</? = <pj,j = 1, -,N} 

These lines are separatrix curve^]: there are no bounded trajectories crossing these boundaries. In the 
field theoretical context, static solutions are also enclosed in these domains. 
Using polar coordinates, the energy functional reads: 

2 



IfdRV 1 R2 (*A' — i 



(12) 



One can think of this functional in (|12j) either as the static energy of the scalar field or, alternatively, 
as the action functional for a particle of unit mass and (R, if) position coordinates moving under the 
influence of a potential —(*U(R,ip) = — (f(R) + jpg((f)) with evolution parameter x. Thus, the motion 
equations of this analogous mechanical system are the ODE system: 

d 2 R = df{R) R2 d^p + dRahp = dg&)_ 

dx 2 dR dx 2 dx dx dip 

Lorentz invariance guarantees that finite action solutions of (|13p are traveling waves of the scalar field 
theory with their center of mass located at the origin; Lorentz and translation transformations provide 
all the solitary wave solutions of the system obtained from the static ones by applying the appropriate 
transformations. 

To solve the ODE (|13p . we take advantage of the existence of two independent first-integrals in the 
mechanical system, namely: 

*-K§) ,+ Hs) , -' w >->> ; • (14> 



technically, these curves are envelopes of separatrix trajectories of the analogous dynamical system. 
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I\ and I2 are respectively the energy and the generalized angular momentum in the analogous mechanical 
system. Written in polar coordinates, the asymptotic conditions (|7|) guaranteeing finite energy in the 
field theory model read: 

lim (R(x), t p(x)) = (R i± , C p j± )eCM , lim ^{x) = , lim ^(x) = . (15) 

x^±oo J x^±oo dx x^ioo ax 

Thus, I\{x = ±00) = and Ii{x = ±00) = 0; because I\ and I2 are invariants of the mechanical system 
they vanish for every point x £ M. Therefore, solitary wave solutions of scalar field theory are in one-to- 
one correspondence (modulo Lorentz boosts) with the trajectories of the mechanical analogous system 
such that I\ = 0, I2 = 0, i.e., trajectories solving the following first-order ODE system: 

g = (-irvmR) , ^(-^v 7 ^) , 0,18=0,1 . (i6) 

We now state the first result concerning the structure of the solitary wave variety in Type II Liouville 
models: 

Proposition I.- There exist kink or solitary wave solutions whose orbits lie on the separatrix straight 
lines. Kink orbits are of two kinds: 1) those connecting two elements (Ri,(fj) and (-R«+i, (fj) of Ai 
through a straight segment (R(x), (fj) (angular rays in the Cartesian internal plane) and 2) those joining 
the vacuum points (Ri,(fj) and (Ri,(fj+i) through the straight line (Ri,(p(x)) (circles on the Cartesian 
internal plane). 

The proof is easy: note that the vanishing of the first integrals I± = 0, 12 = along the kink orbits 
and the assumption of the existence of solutions whose orbit is R = Ri or 99 = (fj on the polar internal 
plane, Ri and cpj being respectively roots of the functions f(R) and g(y>), are compatible conditions. We 
distinguish two possible cases: 



If R — Ri, 



'-5* (£)'->> 



holds and the problem of finding the time-schedule or kink form factor is solved by the quadrature: 



• If if = (fj, the two invariants become 

f ' = K^) 2 - /(fl) = ' ' 2 = ^» = ° ■ 

The vanishing of I2 turns into a identity but the first-order equation annihilating I\ provides the 
kink form factor by means of the quadrature: 

To obtain all the separatrix orbits -those for which I\ = and I2 = 0, in Type II Liouville models- 
it is convenient to apply the Hamilton-Jacobi procedure because the HJ equation is separable in polar 
coordinates and all the trajectories can be found by quadratures. The generalized momenta are pr = ^ 
and p v = R? , whereas the mechanical Hamiltonian (or Hamiltonian energy density in the field theory) 



is: 



n = h R + jgh 9 , h R = -p 2 R - f(R) , hp = ^pl - g{ip) 
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The ansatz of separation of variables J = Jr{R) + Jtp{p) — i\x for the Hamilton principal function 
J converts the PDE Hamilton-Jacobi equation 



into the ODE system: 

2 



1 ( djR \ 

2 V dR J 



dj 
dx 



n 



+ H 



dj_ dj_ 

dR' dip' 



R,ip 



1 f djy 
2\d(p 







9(f) 



-1-2 



k,i2 e 



(17) 



Therefore, the Hamilton characteristic function is given by the quadratures 

, JipW) = sign^) y/2 \ dp^J g{p) - i 2 , 



J R (R) = sign(p R ) V2 / dR\ f(R) + h + 



R 2 



and the kink solutions comply with the equations: 



"Orbit" equation: 



"Time schedule" : #£ 

Oil 



11=0=12 
11=0=12 



7i, 



dR 



R 2 y/V{irj 



sign(p ¥ 



dip 



7i 



(18) 



72, 



sign(pi 



dR 



VvW) 



x + 72 



(19) 



To unveil how the special orbits R = Ri and <p = pj are hidden in the family of trajectories 
parametrized by 71 and 72 , it is useful to show explicitly the connection between the systems of equations 
(|16p and (|18p -(|19 p . To fulfil this goal, we derive from (| 16|) three identities: 



dR 



dR 



x 



{-l)"R? 

dp 



dx 



R?{x) 
dx 

mx] 



+ c 

+ d 



(20) 
(21) 
(22) 



It is clear that (|20|) is equal to (|19|) , whereas subtracting (|22[) from (plj) one obtains (|18|) if 71 = c — d 
and sign^^) = (— l)' 3 . The straight line orbits arise when I71I = 00. In this limit, the (I18p -(|19p system 
makes sense only if (a) p{x) = pj is the orbit, g(pj) = 0, and the time schedule (kink form factor or kink 
profile) is obtained by integrating (|19|) ; (b) R(x) = Ri is the orbit, f{Ri) = 0, and the time schedule is 
obtained by integrating (I22p . Thus, the special straight line orbits are at the boundary of the family of 
kink trajectories. 

Because the kink energy is the action of the associated separatrix trajectory in the analogous me- 
chanical system -| Jr\ + \J v \ computed along the kink path- we state: 

Proposition 2.- The energy associated to a kink or solitary wave solution in the Type II Liouvillc 
models is: 



£[<Kx)] 



dR^/2f(R) 



+ 



d<py/2g(p>) 
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Here, ttr and tt^ are the projectors onto the R- and <p-ax.es in the polar cylinder M + x S ; application 
of these projectors to the x-parametrized kink paths allows us to trade the path integration along com- 
plicated curves by the sum of integrations along straight R(x) = R{ and f>{x) = cpj lines. This is the 
rationale underlying the kink energy sum rules: all the kinks or combinations of kinks having the same 
projections to the polar axes carry the same energy. 

To close this Section on generalities, we briefly describe how the analogous mechanical system is 
related to supersymmetry. In supersymmetric classical mechanics all the interactions are derived from a 
superpotential W(0i,<fo) L28, 29J, related to the mechanical potential energy through the equation: 



1 fdW\ 2 1 fdW^ 2 



+ 2 



(23) 



This is no more than the Hamilton-Jacobi equation for i\ = of a mechanical system with flipped 
potential energy V(0i,<fo) = —U((j)i, fa) (precisely as in the analogous mechanical system), and the 
superpotential is tantamount to the Hamilton characteristic function. The mechanical action reads: 



E 1 r + 

2 J \ dx dx 
and this can be arranged a' la Bogomolny 



dfodcfo dW_dW_ dW dW \ 

dx dx d(f>i d(f>i d(p2 <902 J 



E 



dx 



dx dcf) 



+ 



dfo dW\' 
dx d(f>2 J 



+ 



, dW 



dW 



Thus, the Bogomolny bound Eb = If dW\ is attained by solutions of the first-order differential equations: 

«% _ dW_ dxh_dW_ 
dx d(j)i dx dcp2 



(24) 



Because a first-order ODE system such as (I24D is easier to solve than second-order motion equations, 
it is important to know when the PDE equation fl23f) is solvable. Hamilton-Jacobi separable systems 
admit a complete solution of such an equation and for Type II Liouville systems the situation is: 

Proposition 3. Type II Liouville models admit four superpotentials. 

Proof. In polar coordinates (|23p reads: 



dW 



+ 



1 raw 



dR J ' R 2 
Searching for solutions of (f25|) such that 



dip 



f(R) + ^g(<p) 



d 2 W 







(25) 



(26) 



we plug the expression W(R,tp) 
satisfy the differential equations: 

dF 



_ = (-DV27(S) 



dRdip 

F{R) + G((p) into the previous formula to find that F and G must 

dG 



_ = (-l)V2iM 



with a,/3 = 0,1. Thus, the four superpotentials, the complete solution of the HJ equation for mechanical 
energy equal to zero, are the quadratures: 



W{R, <p) = (-1)" J dR ^/2f{R) + (-if J dip y^icpj 
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The associated first-order equations are our old friends 



dR_dF_ —- dip _ 1 dG ^ yggg 

Note that there is the need of introducing a metric factor in polar coordinates, which we have shown 
to be equivalent to equations (|18 p -(|19 p . obtained through the Hamilton- Jacobi method. It is worthwhile 
mentioning that a global change of sign in the superpotential trades the kink 0k (^0 for the antikink 
solutions 4>k(— x). 

A bonus of the use of the concept of superpotentials is that the two invariants of Type II models can 
be written in Cartesian coordinates in the unified way: 

/ - i f^iV+i f^Y_ i (®l V-- (— V 

1 2 V dx ) 2\dx ) 2\d4i) 2\dfa) 
1 / #i d0 2 \ 2 1 / , dW , dW\ 2 

The vanishing conditions Ji = 0, I2 = 0, required for the kink or solitary wave trajectories, are 
guaranteed by two, rather than one, systems of first-order ODE equations: 



• One expected, equivalent to the ODE system (|24j) : 

d(pi _ dW d(p 2 _ dW 

dx d(pi dx dcp2 

• A new one, unexpected and awkward: 

d4x 0f - 4 dW 20x0 2 dW 



dx 4>{ + 4>\ dcf)i 0f + cf)t, <90 2 

d<\> 2 20102 dW (\>\-(\>\dW 

dx 0f + (fy 001 0f + 0| 502 



^2(01,02) ■ (28) 



System (I28D . however, can be written as the gradient flow equations of a new superpotential W(0): 

dd>i dW dd>2 dW 



dx 50i dx 902 

The reason is that — = 0, and Green's theorem can be applied because the Type II separability 
condition (f2f)|) in Cartesian coordinates reads: 

, , / d 2 W d 2 W \ f± 1 ± 1 . a 2 VF , <9VF 4 „ 
\a0id0i d0 2 a0 2 / a0id0 2 a0 2 a0i 

i.e., precisely the conditions necessary for the curl of the vector field -F(0i, 02) = -^i(0i, 02)ei+-F2(0i, 02)^2 
being zero. 

3 The Al Model 

We shall first choose and 10 W to build the potential energy density: 

c/^)(0i,0 2 ) = (0 2 + 0i-i) 2 + aHl 



\+<t>l) 2 
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Figure 1: Potential energy — U^ A1 \<pi, <fo) in the analogous mechanical system of Model Al: a) Global perspective 
and b) Section showing the singularity and the (f>2 = path. 



depicted in Figure 1. The second summand in the potential energy density U^ifa, fa) spoils the SO (2) 
symmetry preserved by the first summand. The manifold M. of zeroes of C/( A1 )(0i,0 2 ) is a discrete set 
of two elements: 

M = {4+ = (1,0); A. = (-1,0)} . 

The symmetry group of this model is the vierergroup G = Z2 x Z2 generated by the reflections — > —(pi 
and fa — > — (j>2 m the internal plane. This symmetry is spontaneously broken to the Z2 subgroup 
generated by fa —4>2 through the choice of one of the two zeroes, because Ai is the orbit of the 
other Z2 subgroup, in this case generated by (j>\ — * — (j>\. The moduli space of zeroes, the quotient space 
M = M/G = e, however, contains only one element A = {A + ,A-}. 
The partial differential field equations are: 



U 2 



& 2 



dt 2 



dx 2 
d 2 4>i 



1 



dt 2 dx 2 



4-4<j>i 



b 2 + 



a 



2 J.2 



2a 2 



b 2 + 



? 2 ) 2 



+ 



ia 2 ^ 



In the search for solitary waves, however, we follow the general procedure described in the previous 
Section because the analogous mechanical system is a Type II Liouville model with: 



f(R) = (R 2 - 1) 



/(¥>) 



2 ■ 2 
a sin ip 



(29) 



Back in Cartesian coordinates, the superpotential, obtained from the solution of the Hamilton-Jacobi 
equations, reads: 



b 2 + 4> 2 ) 



1 



-l 



J. 



a<Pi 



There are no finite action trajectories leaving the area bounded by the circle G\ = + 4>2 = 1 °f 
unit radius. Also, finite action trajectories do not cross the radial ray r% = <p2 = 0; for this reason we call 
these curves separatrix curved. 



3.1 Solitary waves in the boundary of the kink moduli space 

From Proposition 1 we know that there exist singular solitary waves whose orbits are C\ or R\. We use 
Rajaraman's trial orbit method in order to identify the form factor of these solitary wave solutions. 

2 Do not confuse with separatrix trajectories, those arising for a particular choice of the parameters of the Hamilton 
principal function lying at the frontier between periodic and unbounded motion. Separatrix curves are themselves separatrix 
trajectories in this sense, but very special ones forming the envelop of the whole manifold of these critical motions. 
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• Kl^ 1 " 4 : Plugging the orbit C% = <j)\ + 0| = 1 m t° the field equations the following solitary wave 
solutions are found: 

cj) Klj (x) = tanh V2ax e\ ± sechV2ax e*2 . (30) 

We refer to them as Kl^ using the notation established in Reference [23]; the subscript stands for the 
number of lumps that the solitary wave is composed of - see Figure 2- whereas the superscript specifies 
the elements of Ai that are connected by the kink orbit in a generic way. A can be either A + or and 
by A* we mean the complementary point in M.. In the case depicted in Figure 2, the spatial dependence 
connects the points A + and A- . The density energy and the energy of these kink solutions are respectively 

£ K1 i A *(x) = 2 C j 2 sech 2 v / 2cTX , E(Klf A *) = \W(A-,0) - W(A+,0)\ = 2^2a 

In Figure 2 the main features of the solution are depicted for the choice of plus sign in (|30f) and x = x: 
(a) the kink form factor shows that both field components have non-zero profiles, (b) The kink orbits 
in the internal plane connect the points in M. (c) The energy density is localized around one point and 
we interpret these solutions as basic solitary waves. From a physical point of view they are seen as basic 
traveling particles. 




Figure 2: Solitary waves Kl ± : a) Kink form factor, b) Kink Orbits and c) Energy Density. 



• K2 AA ': Next, the 0i-axis R\ = 4>2 = is tried in the ODE system (|24p . The solution 

(p< 2 i A ( x ) = tanh \/2x ei 

corresponds to solitary waves going from A^ to A + as x varies from — oo to +oo. Note that these 
kinks pass through the origin in the internal plane, a very dangerous point. One could surmise that 
the singularity at this point endows infinity energy to these solutions. We notice, however, that these 
solutions follow the path (f>2 = 0, where the limit of w^\(f)i,(j)2) is zero, and these kinks do not feel the 
singularity. In fact, the kink energy density 



E K2AA ( x ) =2sech 4 V2x 



E(K2 AA * 



W(A^,fo = 0) 



lim W(e, 4> 2 = 0) 



+ 



lim W(e, <j) 2 = 0) - W(A+, (f> 2 = 0) 



4.V2 



is centered around one point in a more concentrated way that the previous solutions - whereas the energy 
itself is different - and K2 AA kinks belong to another class of basic particles in the model. We emphasize 
the following subtle point: because W(4>\, fo) is a regular function all along the K1 A kink orbit Stoke's 
theorem can be applied and E[K if ] depends only on the value of W at the starting and ending points 
A±. The K2 AA * kink orbit, however, hits the origin coming through the <f>2 = axis. The gradient 
of W is undefined in the neighborhood of the origin along this path and careful application of Stoke's 
theorem gives the formula for E[K2 AA *} written above. Note that E[Kl AA *] > E[K2 AA *] if a > § and 
-Ef-fQ^ 1 " 4 *] < E[K2 AA *] if a < |. Figure 3 shows the main features of one of these solitary waves, form 
factor, orbit and energy density. 
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In order to gain a deeper understanding of the physical consequences of passing the kink profile 
through the origin , we shall analyze the stability of K2 AA kinks in some detail. The second-order small 
fluctuation or Hessian operator around these kinks is the diagonal matrix differential operator 



H[K2f A * 



Tin T~(-i2 



. d 2 
dx 2 



4+ 12 sinhivpr 
cosh V2a; 





_ d 2 

dx 2 



-4 + 4 





sinh 2 \/~2x 
cosh 2 \/2x 



I 2(7^ cosh 4 \/2x 
sinh 4 v'22: 



The entry Tin rules the behavior of the tangent perturbations to the kink orbit. Tin iS an ordinary 
Schrodinger operator of Posch- Teller type with Ao = as the lowest eigenvalue. The associated eigen- 
function (or zero mode) describes a perturbation that is a translation of the kink center, see Figure 3(d). 
H22 regulates the orthogonal fluctuations to the 4>2 = orbit. Here, the potential is a positive definite 
infinite barrier with the sing ularity at x = = ^ K2 ^ 4 (0) = 0. This means that fluctuations pushing the 
kink profile away from the origin cost infinite energy, see Figure 3(d). Therefore, the kink profiles of this 
kind of solitary waves behave as strings pinned at the origin of the internal plane. In sum, because the 
spectrum of the Hessian operator is non-negative these solutions are stable. Moreover, Kl^ and ¥^2 AA 
kinks join the same vacuum points asymptotically and in models without this kind of singularity they 
would live in the same topological sector of the configuration space. In the present system, however, they 
belong to different sectors of C because they cannot be homotopically deformed into each other under 
the restriction of finite energy. 



-7.5 -5 -2.5 




2.5 5 7.5 




Figure 3: Solitary waves K2 AA : a) Kink form factor, b) Kink Orbits, c) Energy Density and d) Potential Well 
and Barrier in the two diagonal matrix elements of the Hessian operator. 



3.2 Solitary waves in the bulk of the kink moduli space 

In this model, the Hamilton- Jacobi formulae (jTHJ) and (|19p . or equivalently, the Bogomolny first-order 
ODE system (|271) . become : 



dip 



asm ip 

dR 
R 2 -l 



V2 
V2x 



Integration of (I31D gives the kink orbits: 

(-1)° 

log 



R K (x) - 1 



I - log 


tan y ' 




2 



y/2 



71 



The kink form factor (kink profile) in the R- variable is obtained by integrating (|32|) : 

R K (x) = tanh [V2\x\ 
Plugging this solution into (|31|) . the kink profile for the angular variable is found: 



(31) 
(32) 



(33) 



(34) 
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ip K (x) = 2arctan exp (^/2a(x + 71) — cxcotana; 



(35) 



K AA (j l ,7i~): In Cartesian coordinates, these solutions -made out of two independent pieces- read: 



,K3 AO 



K3 AO , 



tanh Vzlxl tanh 



a 



(V2(x + 



7i J 



cotanh y/2a 



ei 



+ (-l) a tanh\/2|x| sech ct(a/2(x + 71) - cotanh V2x) 



e2 



starting asymptotically from one of the minima denoted as A and and passing through the origin when 
x = 0. This point is reached by each member of this family of solutions in such a way that: 



lim — - — 

s-»0- ax 



lim — — — 

x^o- dx n 







n 



0,1,2,. 



This point is crucial, because every solution of this type can be continuously glued with any solution 
leaving the origin with the same tangency properties and ending at the other point of the vacuum orbit. 
The whole solution on the real line is thus a kink avoiding the singularity, in the same way as the K2 AA 
kink does. In sum, there is a two-parametric family of kinks because we can freely choose HJ trajectories 
with different integration constants 7^ and 7+ in different regions: x < or x > 0. 
A special member of this family is the case when 71 = —00: 



2K^(7i,-oo) 



(x) 



— tanh y/2 \x\ tanh[<r(v / 2(a: + 71) — cotanh y/2x)] e\ + 
+ (-l) Q tanh s/2\x\ sech[cr( v / 2(x + 71) - cotanh s/2x)\ e 2 

tanh y/2 x 



if x < 



if x > 



(36) 

which is formed by gluing a kink in the bulk determined by a finite value of 7 1 at the left of the origin with 
a kink in the boundary at the right of the origin. Note that for this kind of solution the Hamilton-Jacobi 
procedure works independently on the left and right half-lines and this is the reason for the dependence 
on 7]" and 7^. In Figure 4 we depict the kink field profiles for an (7^, —00) orbit with x = and a = 0. 
Several kink orbits are also drawn to show that the peculiarity of choosing 7^ = — 00 changes the ending 
point: these kink orbits link the point A of the vacuum orbit with itself by means of a trajectory crossing 
the origin. 





Figure 4: Solitary wave family K^ lA (7 1 , —00): a) Kink form factor and b) Orbits for the values -f 1 — 0, ±1, ±2. 



Moreover, a Mathematica plot of the energy density for several values of 7j~ (see Figure 5) shows that 
these solutions are indeed a non-linear superposition of the basic solitary waves discussed above: one 
Klf A kink plus one K2 AA kink. Observe, however, that the solutions parametrized by 7f are composed 
of one Kl AA and one-half K2 AA kinks whereas the choice of 7^ = —00 continuously glues the remaining 
one-half K2 AA kink. This justifies the nomenclature K AA (^ , —00) for this kink family. If 7^ is positive, 
these solutions behave as two separate lumps whereas if 7f is negative the two lumps sit on top of each 
other. 
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Figure 5: Energy density of the (71, —00) solitary waves for decreasing values of j 1 



The energy of these solutions is evaluated using expression to find: E[K AA ("f 1 ,—00)] = + 
2y2a. The following kink energy sum rule 



E[Ki A ( 7l -, -00)] = E[K1? A ] + E[K2 AA ] 

holds, showing the composite character of these solitary waves built from two basic kinks. 
Fully general solutions have the form: 



/ j<ao _ _ K AO _ 

^i A \itnT)(^ - ) ^i 3 ( x ! 7i ) ei + <^ 2 3 W7i )e 2 if x < 
l<Ai 3 (a;;7i )ei -02 3 (^;7i )e 2 if x > 

Their energy density is the sum of two Kl^" 4 and one K2 A lumps, and is therefore a composite of three 
basic kinks: 



E[K AA *} = 2E[K1 AA *] + E[K2 AA '} = AV2a + 



4V2 



Figure 6: Several K^ A (j^ ,7j ) kink orbits 

A short remark on stability: K AA (^ , —00) solitary waves are non-top ological whereas the K A (7^ , 7i~)) 
family is formed by topological kinks. Both types, however, are stable because, passing through the ori- 
gin, they cannot decay to lighter kinks with the same asymptotic behavior for identical reasons that 
forbid the decay of the K2 AA to the Kl AA kink. 



4 The Bl Model 

In this section we study the deformation induced on the potential v^ B > (<pi, $2) (HJ) by w^>(<pi, <fo) ©• 
The dynamics is governed by the action functional ([5]), with potential energy density: 

U^\4>iA2) = (0? + 4~ l) 2 (^ + <Pl ~ a 2 ? + . 

{<Pl + 92) 

Thus, as in the Al model, by adding the if/ 1 ) perturbation the symmetry under the SO (2) group is 
explicitly broken to the discrete subgroup G = Z2 x Z2 generated by reflections of the fields: 4>X ~^ —'Pi-, 

<t>2 -02- 

For the sake of simplicity we set the value of a to be: a = 2. Other non-null values of this parameter 
yield the same qualitative behavior of the system with only analytic differences in the specific expressions. 
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The set of zeroes ) is a continuous manifold with two connected components: the disjoint 

union of two circles of radius 1 and 2, S^ =1 U S^ =2 - The w 1 perturbation, however, forces the set of zeroes 
of t/( B1 ) to be the discrete set of four points: 



M = {B~ 



-2,0) ; A- = (-1,0) ; A+ = (1,0) ; B+ = (2,0)} 



The vacuum orbit is the union of two G-orbits: A4 = {A , A + } U {B , B + }. The moduli space of vacua 
M = M/G = A U B, however, is the union of only two elements: A = {A~,A + }, B = {B~,B + }. 

The field equations are: 



d 2 



dt 2 



d 2 c 



dx 2 



J2±2 



(1 - cj>\ - 0|)(4 - 4>\ ~ $)(5 - 2<p 2 - 2<g) + 



(l-^-0|)(4-0?-^)(5-2^ 



2^) 



+ 



(7 



a 



2x2 



2(0? 



9t 2 dx 2 

The analogous mechanical system is the Type II Liouville system for which: 
f(R) = (l-R 2 ) 2 (i-R 2 ) 2 , g{tp) 



2\3 



2 ■ 2 

a sin <p 



Back in Cartesian coordinates, the superpotentials, obtained from the solution of the Hamilton-Jacobi 
equation, read: 



W{<f> 1 ,(f>2] = V2- 



+ ■ 



1 



l + a 2 



% + 4>l) + a 2 



acpx 



There are no finite action trajectories trespassing the circle: C2 = 4>t + 4>2 = 4- Also, C\ = 4>f + 4>2 = 1, 
and i?i = 02 = are not crossed by any finite-action trajectories. Thus, C%, C2, and R\ form the set of 
separatrix curves of the Bl model. 



4.1 Solitary waves at the boundary of the moduli space 

From Proposition 1 in section §2, the existence of singular solitary waves with C2, C\, and R\ orbits 
follows. Use of the Rajaramans trial orbit method is effective. 

• Kl±* . Substituting (j) 2 + cf) 2 , = 1 into (|27p we obtain the following solitary wave solutions: 

Kl1 (x) = tanh \[2ax e\ ± sechv^crx^ , 

joining the points A + and A^, see Figure 7. They are basic lumps or particles of the Bl model - recall 
that exactly the same solutions are also solutions of model Al- and their energy density is concentrated 
around a point in the real line, see Figure 7c. Specifically, we find that the energy density and the energy 
of these kinks is exactly the same as in model Al: 



£ K1 ? A \x) = 2a 2 sech 2 V2ax , E[Klf A *] = \W(A-,0) - W(A+,0)\ = 2^2a 



• K^ B * . Analogously, we choose (f) 2 + (j) 2 = 4 as the trial orbit and plug this expression into (|27p , to 
obtain: 



K i (x) = 2 tanh — — e 1 + 2 sech 



crx _ „ , ax _ 

e 2 



2y/2 ' 2V2 



15 



Figure 7: Solitary waves Klf A : a) Factor Form, b) Orbit and c) Energy Density. 



These solitary waves share similar features with the K AA * kinks, although in this case the asymptotically 
linked points are B + and Even though the total amount of energy carried out by the K BB * and 
K1 AA * kinks is the same, the energy density is less concentrated in the Kf B kinks: 



„K-i BB* , . 

£ K1 i (x) 



a 2 ax 
-sech 



£[Klf B *] = \W(B-,0) - W(B+,0)\ = 2V2a 



2 " 2V2 

see Figure 8. They describe different basic particles, living in different topological sectors. 
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Figure 8: Solitary waves Kf B : a) Factor Form, b) Orbit and c) Energy Density. 



On the fa = axis there are two different kinds of singular solitary waves. 

• K AB : The first kind includes two kink orbits, the closed intervals -K = 1- and [^4+,S + ] 

= 0-. The corresponding kinks connect points in different elements A and B of the vacuum moduli 

space. We denote this kind of solitary wave solutions as K AB kinks and their profiles or form factors are 

easily found to be: 



) K i (2) = 2 cos 



<o\f2x 



2vr. 



- arctan e uv ^ H K 

3 3 



ei 



« = 0,1 



The energy density and total energy of these solutions is: 

£ K i B (x) = 32sech 2 6v / 2x sin 2 Q arctan e 6v ^ s + ^-i^j 



E{K AB ] = \W(A T ,0)-W(B T ,0)\ 



22^/2 
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The energy density is again localized around a point and these solutions are also basic solitary waves or 
lumps of the Bl model, see Figure 9. 

• K2 AA * : The orbit of the second kind of solitary waves on the 4>2 axis is the L4_,A + ] interval. The 
kink profile of these kinks connects these points in the vacuum orbit and is also easy to determine: 



<A K2 ^*(x) = 2cos 



6V2x 



47T 



- arctan e ov zx + — 



ei 
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Figure 9: Solitary waves ; a) Form Factor b) Orbit and c) Energy Density. 

Despite being analytically different, the behavior of these kinks is completely analogous to that of the 
bizarre K2± solitary waves of the Al model. The energy density and total energy of these kinks are: 

S K2 * A * (x) = 32 sech 2 6\/2x sin 2 Q arctan e 6 ^ + y ^ 



E[K2{ 



AA*i 



W(A_ ,^2 = 0)- lim W(e, <t> 2 = 0) 

£-+0" 



+ 



lim W(e, 4> 2 = 0) - W(A+, (f> 2 = 0) 



76V2 
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see Figure 3 to find qualitative plots of their properties. 



4.2 Solitary waves in the bulk of the moduli space 

The Hamilton- Jacob i orbit (|18p and time schedule (|19p are in this case the quadratures: 



dR 



R 2 ){A-R 2 ) 



(-1)/ 



cr sm 99 



V2x 



{l-R 2 ){A-R?) 

Integration of these equations provides the analytic expressions for kink orbits and form factors: 



(-l) f 



4^ 6 (R + 2)\(R-\) 



-I) 13 <p 

log tan - 

a 1 



12 S (R + 2)(R- l) 2 



V271 
\/2x 



(37) 
(38) 



Solitary wave solutions only arise from equations (|37h and (138f) ifl<i?<2orO<-R<l. Let us 
denote G = arctan e 6v/2x + in the first range, G = arctan e Gv ^ x + 4^ in the second range, and let us 
define the functions: 



fii(x) 



cos f + \ 
cos^-i 



giCg) - 1 1 2Q 
sMx) + 1 2 3 



A(x) = tan 



G 



(39) 



• K^ B (7i): If 1 < R < 2, the solutions in Cartesian coordinates are: 

^ B M(x) = ^ B ^\x)e l+ ^ B ^\ X ) 



(-2 



x) 



1 



n 2 (x) 
+ (-IT 



l + e -2V2a 7l e -f na(x) n 8 A ^ ^) 
2 e V2a7l e |n2(s) Of(x)A^(S) 



ei 



n 2 (s) 



^ + n 1 s (x) Ah(i) 



e 2 



a = 0, 1 
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The kink orbits belonging to these two one-parametric families connect either A_ to B + if N = 0, or A + 
to B- if H = 0, and are confined inside the annulus bounded by the circles C\ and C 2 , see Figure 10. The 
energy density of each member of these kink families is localized around two distinct points, see Figure 
10. For sufficiently negative, values of 71 the two lumps are composed of one K1 AA * and one K AB basic 
kinks. In the opposite regime, for sufficiently positive values of 71, two lumps arise again, close to one 
K BB * and one K AB basic kinks. Note that the basic kinks appear at the boundary of the moduli space 
when either 71 = —00 or 71 = 00. Intermediate tuning of 71 continuously shifts from one configuration 
to the other, see Figure 11. These features are synthesized in the notation ¥^ AB (^\). 





Figure 10: Solitary wave family (71)-' a) form factor and b) Orbits for values 71 = 0, ±1, ±2. 
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Figure 11: Energy Density of the solitary wave family K A (71) for increasing values 0/71. 

• K^ A (7 1 f ,7jf): Inside the C±, < R < 1, there are kink solutions completely analogous to the 
^■3 A (li i li) solitary waves of the Al model, see Figure 6. The kink profiles have the form: 



^•(7+71-) 



(x) 



■ (s; 7 X 



,K AO 



K AO _ 

e*i + 4> 2 3 (x; 7 X ) e 2 if x < 

K AO + 
(x;7i) e*i - 4> 2 3 (x; 7^) e 2 if x > 



where 



^°(7l) 



(X) = 



1 



n 2 (x) 
(-i) Q 



- 1 



fi 2 (x) (e 2 v^7i e f + (x) An(i)J 



ei 



e 2 



Q: 



0,1 



and 6 = arctan e 



6a/2S 



+ ^ in the < R < 1 range. 



In particular , the limiting case K AA {^)^ , —00) gives a solitary wave identical to the K AA (^fi) kink of 



K AO _ K AO _ -i„„ M . 

the Al model. To obtain the kink profile one replaces ^ 1 3 (x; r y l ) ei+0 2 3 (x; r y 1 ) e 2 by <f> Kl ^ (7i)(x) = 
2cos6e*i on the left (x < 0). We recall that these solutions connect a point in A with itself through a 
trajectory that crosses the origin, see Figure 4. 

In sum, the structure of the solitary wave variety in the Bl model is as follows: There exist four 
kinds of basic lumps associated with the Kl^ 1 " 4 , ¥±1 AA , Kf B and K AB kinks, which display different 
distributions of one-point localized energy densities. The K^( 7 +, 7l )), K^( 7l ), and K^ B ( 7l ) kinks, 
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however, are composite solitary waves. The K^ A ( 7l + ,7 1 ~) kinks are combinations of two Kl^" 4 and one 
K2 AA basic lumps; the K^ j4 (7i) are formed by one K1 AA * and one K2 AA * basic lumps, the K^ B ( 7l ) 
kinks exhibit an orbit dependent structure that varies between the limiting combination formed by either 
the Klf" 4 and K AB kinks or the Klf B and K AB kinks. Therefore, the following kink energy sum rules 
hold: 

E[K£ A } = 2E[K1 AA ] + E[K2 AA ] = A^2a + — V2 

E[K AA ] = E[K1 AA ] + E[K2 AA ] = 2^2a + — V2 

15 

E\K AB ] = E\K AB ] + E[K BB ] = E[K AB ] + E[K1 AA ] = 2^2a + —V2 . 

15 

With respect to stability, all the kinks inside the C\ circle behave like their cousins in the Al model. All 
the kinks in the 1 < R < 2 annulus are stable. 



5 The A2 model 



In this Section we analyze the deformation induced on the potential v^ A \4>i, ^2, 
The action functional ([5]) with potential energy 



© by w 



(2)/ 



U A 2 



l) 



1) 2 + 



a 



(40) 



governs the dynamics. This deformation explicitly breaks the SO(2) symmetry to the discrete group 7Li 
generated by the reflection <p2 —> —<p2- The set M of zeroes of U((f>i, 4>2) has only one element: 

M = {A = (1,0)} 

In this model, the second order partial differential equations 



2 



2 



dt 2 



dx 2 
d 2 6. 



dt 2 dx 2 



+ 



a 2 ^-202 + 2^ lv VF + 



+ 



a 



20 2 



30102 



b\ + <t> 2 ) 2 



are the field equations. The asymptotic conditions (J7|) guaranteeing finite energy compel the solitary 
wave solutions to link the only minimum A with itself. The general results in Section §. 2 ensure that 



the Ci = <j){ + 



1 circle and the n = 02 = U 0i > abscissa half axis are separatrix curves. Unlike 



the Al and Bl models (see Section §. 3 and §. 4) the singularity at the origin is of different nature in 
this case. We know from previous Sections that there exists a path 02 = along which the singularity 
is not felt in either model Al or Bl. The potential energy C/ A2 (0i, 02), however, reduces on the abscissa 
axis to: 



^ A2 (0!,O) = (0^-1) 



sign (0i )) 



and the singularity is always felt in the 0i < negative half-axis. There is no escape to kink orbits that 
enter through the positive abscissa half-axis, and finite energy solitary waves are pushed away from the 
origin. 

Nevertheless, there are kink orbits on the C\ = 0? + 0n = 1 circle. 



19 



: Plugging this curve into (I27D we obtain 



rw AA / \ i - i 9 ff% . \ -, „ , ax , ax _ 
(fi^i (x) = [2 tantr — - 1 ei + 2 sech — tanh — e 2 

v2 / v2 v2 



a solitary wave connecting the point A with itself when x varies from — oo to oo, see Figure 12(a,b). The 
energy density 

E[K AA ] = 4V2a 



£ K ? A (x) = 2a 2 secb 2 ° X 



V2 



is localized at one point so that the K AA is a basic lump, see Figure 12(c). 
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Figure 12: Solitary waves K AA : a) Factor Form, b) Orbit and c) Energy Density. 

Again we are dealing with a mechanical analogous system of Liouville Type II such that: 
f(R) = (i? 2 - I) 2 , 5 M = a 2 sin 2 | . 

In Cartesian coordinates the superpotentials read: 



W+(fa, fa) = (fijtf - l) - 2a(-Wl + y^ 



W 



(-i r V2y^T^ 2 ( ^ - 1) + 2a{-l)\ l + 



and we realize that the kink energy is a topological bound 

E[K AA ] = \ W + (A) - W~(A) \ = AV2a 

giving the winding number of the kink orbit around the origin. This fact ensures stability for the K AA 
kink. 

The Hamilton- Jacobi theory yields the remaining solutions confined inside the unit circle C\. In 
this model, the Hamilton-Jacobi formulae (|18p and (|19p or, equivalently, the Bogomolny first-order ODE 
system (|27|) become: 



-1) 



dR 



R 2 {R 2 - 1) 



dtp 



2y/a sin ^ 



dR 



R 2 -l 



V2 7l 
V2x 



(41) 
(42) 



Integration of (|4ip gives the kink orbits: 

, / R K (x) + l 
log 1 



R K (x) - 1 



2 1 



| - log 


tan^ ; ; 
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a/2- 



71 



(43) 



The kink form factor (kink profile) in the R- variable is obtained by integrating f|42|> : 

R K (x) = tanh [y/2\x\ 
Plugging this solution into (|4ip the kink profile for the angular variable is found: 

(f K (x) = 4arctan exp (V2a(x + 71) — ercotanV^a 
Back in Cartesian coordinates, we obtain: 



tanh 



\/2>| { 



1 — 2 sech 



2 (71 



V2x + cotanh VZx^J j ei 



+ 



(44) 



(45) 



+ 2 tanh v2| x | tanh ( 71 — V2x + cotanh y/2x ) sech ( 71 — \Fix + cotanh \pix ) e*2 



Some orbits for these solutions are depicted in Figure 13. Here, we notice that all solutions of this kind 
head towards the origin; thus, they cannot be regarded as finite-energy solutions. 



Figure 13: Trajectories of solutions confined in the unit circle C\. 

In sum, in model ^42 there exists only one basic solitary wave, the Kf~ A kink, and no composite 
solitary waves arise. 

6 The B2 Model 

Finally, in this last section we study the deformation u/ 2 ) ([9j) induced on v^ B ' (jlj). Therefore, the potential 
term 

v B 2(<h,<h) = (<Pi + 4- 1) 2 (0? + 4- a 2 ) 2 + ~ ( j\,2, ( 1 " -7=f=f ) (46) 

enters the functional action ([5j). We set the a = 2 value bearing in mind that models of B2 type 
characterized by other non-null values of the parameter a show a very similar behavior. As in the A2 
model, the second summand on the right hand side of (|46f) explicitly breaks the SO (2) symmetry of 
the unperturbed system to the discrete group Z2 generated by the 02 — > — 02 reflection. The non- 
deformed manifold of zeroes M. a =o = Sj. =1 U S^ =2 , the disjoint union of two circles, becomes the discrete 
set of two elements: 

M = {A = (1,0) ; 5 = (2,0)} 
The Euler-Lagrange equations form the system of second-order PDEs: 



- -4>i~ 2 )(4 - <Pl - P2){5 - 2( t>2 ~ 2( t>2) + 



dfi dx* — - — - — "™ ' 2 (0 2 + 0|)§ 



202 30102 



b\ + 02)2 (02 + ^ 
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Solitary waves are solutions of these equations complying with the asymptotic conditions (J7|). Thus, the 
kink must connect a point, either A or B, with itself or must link A with B. The general results in section 
§. 2 ensure that the separatrix curves in this model will be the C\ = (f>\ + </>f = 1> C*2 — + 0| = ^ 
circles, and the r\ = <p2 = U <p\ > abscissa half-axis. The nature of the singularity arising in the 
(|46|) potential is identical to the singularity of the A2 model, see section §. 5. Each solution that crosses 
through the origin cannot be regarded as having finite energy. 

The analogous mechanical system is a Liouville Type II system such that: 



f(R) = (R 2 - lf{R 2 -Af 



( \ 2 ■ 2 V? 

g(cp) = a sm - 



Accordingly, in Cartesian coordinates the superpotentials read: 



1(0? + 01) + 4] 



w-iyUi + 



62)=V2(-l)«vW+^[l(^i + 



<2\2 



^? + 0i) + 4] +2^-1)^/1 + 



<t>\ 



<i>\ 



h > o 
h < o 



Bearing this in mind we now describe the solitary wave or kink variety of the model: 
6.1 Solitary waves at the boundary of the moduli space 

• K±- A : Proposition 1 in section §. 2 guarantees the existence of solitary waves with orbits running on 
the unit circle C\ = (f>\ + <^>| = 1. Integration of the first-order equations (I27p for configurations living on 
C\ provides the solitary wave solutions: 



2 cx \ ^ ax ax . 

2 tanh —= — 1 ) e± + 2 sech —= tanh —= e2 



V2 J V2 ' V2 

These solitary waves are non-topological -but stable- kinks identical to the K AA kinks of the A2 model, 
encircling the origin in a orbit that starts and ends at A, see Figure 14(a,b). They are of course basic 
lumps because their energy density 



2 a 2 sech 2 ^ 



V2 



E[K AA ] = \ W + {A) - W~(A)\ = 4V2a 



is localized around a point, see Figure 14(c), and their total energy is a topological bound proportional 
to the winding number of the kink orbit around the origin. 





1.5 

1 



Figure 14: Solitary waves : a) Form factor, b) Orbit, and c) Energy Density. 



• Kf B : By the same token there are solitary waves whose orbit is the C2 = <fi± + 4>\ = 4 circle. 
Plugging this trial orbit into the first-order equations we obtian via one quadrature the solitary wave 
solutions: 



^? B (x) = (4 tanh 2 



ax \ ^ ax ax _ 

— -= — 2 ei + 4 sech — — tanh — — eo 
4^/2 / 4^/2 4V2 
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These non-topological kinks have orbits departing from and arriving at the same point B, see Figure 
15(a,b). Their energy density 



S K? B ( X ) = iLsech 2 -1±= , E[K BB ] = \W + {B) - W~(B)\ = 4y/2a 

2 4\/2 

is again localized at a point, but the energy distribution is different from the energy distribution of the 
K AA kinks. Thus, the K BB kinks, besides belonging to a different topological sector of the configuration 
space, are basic lumps or particles different from K AA kinks, see Figure 15(c). Moreover, K BB kinks are 
also stable because their energy is given by the winding number of the kink orbit around the origin. 




Figure 15: Solitary waves Kf B : a) Form Factor, b) Orbit, and c) Energy Density. 



• K AB : In this model there also exist kink orbits on the separatrix curve r± = fa = U fa > 0. 
These kink orbits connect the points A and B. Plugging the abscissa half-axis as a trial orbit into the 
first-order equations (|27p . the following solitary waves are found: 



TK AB / | 2 

4> 1 (x) = 2 cos I — arctan e 

These kinks correspond to a third kind of basic lump or particle in the B2 model. Their energy density 

22^2 



£ AB (x) = 32sech 2 6y/2x sin 2 Q arctan e^^j 



EK 



ABi 



\W ± (E) - W ± (B)\ 
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is concentrated around a point in a different fashion to the energy densities of the iso-energetic K AA and 
K BB kinks, see Figure 16. E[K AB is an absolute minimum of the energy in the C AB topological sector 
and these kinks are stable. 



i 

0.5 




Figure 16: Solitary waves : a) Form Factor, b) Orbit, and c) Energy Density. 



6.2 Solitary waves in the bulk of the moduli space 

• K AB (71 ) : Finally, by applying the Hamilton- Jacobi procedure we find the generic solitary wave solutions 
in the bulk of the kink moduli space. There are solutions confined inside the unit circle C\ that go to 
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the origin. They are identical to the infinite energy solutions of the A2 model and must be discarded as 
solitary waves. 

There are, however, genuine solitary wave solutions confined inside the annulus delimited by the C\ 
and C2 circles. These kink orbits connect the two points in M, A and B. In this case ,the Hamilton- Jacobi 
orbits (fTHl) and time schedules (fl~9l) are provided by the quadratures: 



(-1) C 



dR 



R 2 (l- R 2 )(4- R 2 ) 



-1) 



dip 



asm % 



dR 



V2x 



{l-R 2 ){A-R 2 ) 

Integration of these equations provides the following analytic expressions for kink orbits and form factors: 



-1 1, \R-2\s(R+l) 
log 



4R 6 B (fl + 2 )l|fl-l| 



2(-l) /3 1 V? 

log tan -f 

a 4 



1 , \R-2\(R+ l) 2 



V271 
V2x 



12™° (R + 2)(R-1) 2 

if 1 < R < 2. If < R < 1, the orbits above pass through the origin and have infinite energy. 
In Cartesian coordinates, these solitary wave solutions are: 



jj>"<*»>(*) 



^(71), 



n 2 (x) 



1 + 



g e 2V2<T7i e fn 2 (x) 



1 + e-^i e -f ^(x) Q6 ^ A ^ ^ 



e V2<T7i e f a 2(s) + o|(s)AM(x)l 
where we have made use of the notation defined in ()39|) . see also Figure 17. 





Figure 17: Solitary wave family Kj (71): a,) Form factor, and b) Orbits for several values 0/71. 

The distribution of the energy density reveals that solitary waves of this type are composite. For 
sufficiently positive values of the orbit parameter 71 we find a combination of one Kl AB and one K AA 
lump. This structure turns into a combination of one K1 BB and one K AB lumps for sufficiently negative 
values of 71, see Figure 18. For intermediate values of 71 the two basic kinks are entangled. In any case, 
for every member of this family the following kink energy rule holds: 

_ 99 _ 

E[K£ B M] = E[K AA ] + E[K AB ] = E[K BB ] + E[K AB ] = 4^ + £ V2 . 

15 
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Figure 18: Energy density of K2 B (ji) solitary waves for decreasing values 0/71. 
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